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Abstract 

We study the propagation of gravitons within 5-D supersymmetric braneworld mod- 
els with a bulk scalar field. The setup considered here consists of a 5-D bulk spacetime 
bounded by two 4-D branes localized at the fixed points of an orbifold. There 

is a scalar field in the bulk which, provided a superpotential determines the 

warped geometry of the 5-D spacetime. This type of scenario is common in string 
theory, where the bulk scalar field is related to the volume of small compact extra 
dimensions. We show that, after the moduli are stabilized by supersymmetry breaking 
terms localized on the branes, the only relevant degrees of freedom in the bulk consist 
of a 5-D massive spectrum of gravitons. Then we analyze the gravitational interaction 
between massive bodies localized at the positive tension brane mediated by these bulk 
gravitons. It is shown that the Newtonian potential describing this interaction picks 
up a non-trivial contribution at short distances that depends on the shape of the su- 
perpotential We compute this contribution for dilatonic braneworld scenarios 
oc e"**^ (where a is a constant) and discuss the particular case of 5-D Heterotic 
M-theory: It is argued that a specific footprint at micron scales could be observable in 
the near future. 



1 Introduction 



Recent tests of gravity at short distances [H El [HI H] have confirmed that Newton's 
inverse- square law holds down to a length scale 56/im. This has substantially improved 
previous constraints on exotic interactions mediated by the exchange of massive scalars 
or vectors between neutral atoms [5], where a Yukawa type contribution to the New- 
tonian potential is generally expected. It has also lowered the scale at which large 
extra dimensions [HI El E] and braneworld models [HI [IHl [H] may show up by affecting 
the propagation of gravitons in the presence of a large -or infinite- extra dimensional 
volume. 

Indeed, in theories where matter fields confine to a 4-D brane and gravity is the only 
massless field able to propagate along the extra dimensional volume [121 lEl [E] , one 
generally expects short distance corrections to the usual 4-D Newtonian potential. The 
shape and distance at which these corrections become relevant generally depend on the 
geometry and size of the extra dimensional volume, thus allowing for distinctive signals 
dependent of the particular content of the theory. In the single-brane Randall- Sundrum 
scenario [TT] , for instance, where a 4-D brane of constant tension (x k is immersed in an 
infinitely large five- dimensional AdS volume, a zero mode graviton g^i, localizes about 
the brane. This zero mode is exponentially suppressed away from the brane with a warp 
factor oc e~'^^, where z is the distance from the brane along the fifth extra-dimensional 
direction. The Newtonian potential describing the gravitational interaction between 
two bodies of masses mi and m2 localized at the brane, and separated by a distance 
r, is then found to be [IH [IS HEl [HI [18] 



where Gn is Newton's constant. The correction 2/3A;^r^ springs out directly from 
the way in which bulk gravitons propagate in an AdS five-dimensional spacetime. A 
correction like this provides an important signature for the low energy phenomenology 
of braneworld models with warped extra-dimensions; if the tension k is small enough 
as compared to the Planck mass Mpi = (SttGn)"^''^; then it would be possible to 
distinguish this type of scenario from other extra-dimensional models in up-coming 
short distance tests of gravity (present tests give the robust constraint 1/k < 11/im). 

It is therefore sensible to ask how other braneworld scenarios may differ from the 
Randall- Sundrum case at short distances, especially within the context of more realistic 
models. The purpose of this paper is to shed light towards this direction. Here we 
study the propagation of gravitons within 5-D braneworld models where the geometry 
of the extra-dimensional space differs from the usual AdS profile. We will show that 
the gravitational interaction at short distances is sensitive to the geometry of the 
extra- dimensional bulk in such a way that the Newtonian potential picks up a non- 
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trivial correction at scales comparable to the tension of the brane. As we shall see, 
this correction may differ dramatically from the one depicted in Eq. ([T]). We refer to 
[ini EHl ED [221 [231 [211 [2S1 [2S| for other works on short distance modifications to general 
relativity within the braneworld paradigm. 



1.1 General idea 

We will look into a fairly general class of supersymmetric braneworld scenarios with a 
bulk scalar field (f). The model considered here consists of a 5-D bulk spacetime bounded 
by two 4-D branes localized at the fixed points of an S'^/Z2 orbifold. The tensions of 
the branes are proportional to the superpotential of the theory, allowing for BPS 

configurations in which half of the bulk supersymmetry is broken on the branes [2Zl [2H] • 
These types of models are well motivated from string theory, particularly within the 
heterotic M-theory approach [30] where, in the 5-D effective low energy theory, the 
scalar field is related to the size of the volume of small extra-dimensions compactified 
on a Calabi-Yau 3- fold [3ll [32] . To gain insight into the gravitational phenomenology 
of this model, we shall only consider the bosonic sector of the theory. 

One crucial property for us coming from this class of models is that the warping of 
the extra-dimensional volume depends on the form of the superpotential To be 

more precise, given a metric g^u = uj'^{z)rj^y, where ri^^ is the usual Minkowski metric 
and z is the coordinate parameterizing the proper distance along the extra- dimension, 
then bulk fields are related to by means of the following first order differential 

equations 

— = - W and 0' = — -, 2 

where ' = dz- In the particular case of Randall- Sundrum branes, the superpotential 
is simply a constant W = 4k, implying an AdS5 spacetime with u = e~^^. For 
more general superpotentials, however, the warping of the extra-dimension may have a 
richer structure and even contain singularities [33]. For example, in the case of dilatonic 
braneworlds, where W{(j)) = Ae°"^, one encounters a singularity = at a distance 
z = l/a'^W{(f)i) from the positive tension brane, where 0i is the value that the scalar 
field acquires on the brane. In this way, while in the Randall- Sundrum model there is 
a positive-tension brane in an infinite volume (with a; = at infinity), in more general 
cases one may have a single positive-tension brane immersed in a finite bulk- volume at 
a certain distance from the u = singularity. 

The basic idea of this paper is to compute the effect of such geometries on the gravi- 
tational interaction between massive bodies localized on the same brane. For simplicity, 
we shall consider single-brane configurations in which the visible brane corresponds to 
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the positive tension brane, whereas the negative tension brane locahzes at the bulk 
singularity. 

1.2 The moduli problem 

Models with extra-dimensions generically predict the existence of massless degrees of 
freedom, the moduli, at the 4-D effective theory level [IHl [3l] . These moduli appear 
coupled to the matter sector with the same strength as gravity, leading to significant 
long range modifications to general relativity, well constrained by both Solar system 
[35l[36l|37] and binary pulsar tests [3HI EHl HOl HI] . In the Randall-Sundrum model, the 
radion moduli vanishes as the negative tension brane disappears at infinity -the single 
brane limit- leaving gravity as the only relevant long range interaction of the model. 
However, in more general braneworld scenarios -as the one we consider here- one 
typically expects other massless degrees of freedom, even in the single brane limit. For 
instance, current Solar system tests imply a constraint on the following dimensionless 
parameter a = W'^d^pW 

< 1.5 X 10"^ (3) 

where a is evaluated on the brane in which tests are performed [12]. The Randall- 
Sundrum model corresponds to the trivial case = 0, thus passing the test with 
flying colors, nevertheless, in more realistic models one has ~ 0{1). As we shall 
learn later in more detail, in order to have significant effects at short distances -say, 
at micron scales- different from the Randall-Sundrum case, it is necessary to be in the 
range ^ This strongly contrasts with the bound of Eq. (I3l). 

One way out of this problem consists in taking into account a stabilization mechanism 
for the moduli [ISl HI], in this case, for the bulk scalar field (p. If becomes massive on 
the branes, then the only relevant long range interaction in the bulk would consists of 
the gravitational field. To this extent, we consider supersymmetry breaking potentials 
localized at the orbifold fixed points. We will show that, provided certain simple 
conditions on these potentials, it is possible to stabilize the bulk scalar field without 
spoiling the vacuum geometry of the extra-dimensional space dictated by Eq. ([2]), and 
therefore retaining all the interesting features coming from the bulk-curvature. 

1.3 Plan of the paper 

This work is organized in the following way: We start in Section [2] by introducing 
braneworld models with a bulk scalar field. There, we deduce the equations of motion of 
the system and present the zero mode background solution -a BPS vacuum state- and 
its effective theory, a bi-scalar-tensor theory of gravity. The scalar degrees of freedom 
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of this theory consist of the boundary values of (f) at both branesE Then in Section [3] we 
study the hnear perturbations of the fields around the zeroth-order solution presented 
in Section [21 There, we also consider the problem of stabilizing the moduli. We 
show that, once the zeroth-order moduli are stabilized, the only relevant degrees of 
freedom at low energies in the bulk consist of a massive tower of 5-D gravitons. In 
Section m we analyze the effects of these massive states on the gravitational interaction 
between massive bodies localized at the same brane. There, we compute the short 
distance modifications to Newton's inverse-square law parameterized by a function 
/(r) appearing in the form 

V{r) = -Gn^ [1 + /(r)] . (4) 

As we shall see, the shape of the function /(r) depends heavily on the shape of the 
superpotential W{(f)). To put things into context, we further compute the corrections 
arising in the particular case of dilatonic braneworlds W{(f)) oc e°"^, where a is a con- 
stant. There, we also discuss the particular case of 5-D Heterotic M-Theory and the 
prospects of observing these corrections in the near future. Finally, in Section E] we 
provide some concluding remarks. 



2 Braneworld models with a bulk scalar field 

Let us consider a 5-D spacetime M with topology M = x 3^/7^2, where is a 
fixed 4-D Lorentzian manifold without boundaries and S^/7j2 is the orbifold constructed 
from the one-dimensional circle with points identified through a Z2-symmetry. M is 
bounded by two 3-branes located at the fixed points of j^L^- We denote the brane 
hyper-surfaces by Ei and S2 respectively, and call the space M bounded by the branes, 
the bulk space. In this model there is a bulk scalar field with a bulk potential 
and boundary values 0^ and 0^ at the branes. Additionally, the branes have tensions 
Ai and A2 which are given functions of 0^ and 0^, respectively (see Fig. [1]). The total 
action of the system is given by 

'S'tot = S'bulk + 'S'brane, (5) 

where S'bulk is the term describing the gravitational physics at the bulk (including the 
bulk scalar field) 

^bulk = ^ / - ^ / {m + f/(0)] + ^GH. (6) 

* Alternatively, one may define the two scalar degrees of freedom as the distance between the branes 
(the radion) plus just one boundary value of </) at a given brane. 
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Figure 1: Schematic representation of the 5-D brane configuration. In the bulk there is a scalar field 
(/) with a bulk potential U{(j)). Additionally, the bulk space is bounded by branes Ei and E2 located 
at the orbifold fixed points. The branes are characterized by tensions Ai and A2, and contain matter 
fields and ^'2 respectively. 



Here the integral symbol J is the short notation for J d^Xyf—g^, where X^, with 
A = 1, ■ ■ ■ , 5, is the coordinate system covering M and is the determinant of the 
5-D metric gAB of signature ( — h + + +)• M5 is the 5-D fundamental mass scale and 
R^^^ is the 5-D Ricci scalar. Observe that in the present notation the bulk scalar field 
(J) is dimensionless. The third term of Eq. (El) corresponds to the Gibbons- Hawking 
boundary term Sgh oc Jj._^ K — J^^ K, added to make the bulk gravitational physics 
regular near the fixed points. The action ^brane appearing in Eq. ^ stands for the 
fields at the fixed points. It is given by 

^brane = ^matter - ^ / Ai(0l)-^ / \,{ct>'), (7) 

where Ai(0^) and A2(0^) are the brane tensions and ^matter the action describing the 
matter content of the branes, which we write 

^matter = ^7^1^] + ^21^2, 9^}], (8) 

where \E'i and \E'2 denote the matter fields at each brane, and g^^^u and g^l^} are the in- 
duced metrics at Ei and S2 respectively. In what follows we summarize some important 
properties of this system. 

2.1 5-D supergravity 

As already mentioned, we focus our interest on a class of models embedded in super- 
gravity, where the bulk potential and the brane tensions Ai(0^) and A2(0^) satisfy 
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a special relation so as to preserve half of the local supersymmetry near the branes |27j . 
The relation turns out to be 



and As = ~W{(f)'^), 



(9) 



where W = W{(j)) is the superpotential of the system. Observe that the tensions Ai 
and A2 depend on W with opposite signs. 

Several aspects of this class of models have been thoroughly investigated over the last 
few years, among them: Braneworld inflation [IS, HG], HT] , their low energy dynamics 
[m mi Eni EU [52], Uni ISl], brane collisions |55l [56], and various phenomenological 
aspects ISZl EHl EH]- This class of model is attractive for several reasons: On the one 
hand, they offer a natural extension to the much studied Randall-Sundrum model, 
where a fine tuning condition between the bulk cosmological constant and the tensions 
allows a null effective cosmological constant on the brane. This is also the case here 
[601 [611 [621 [63] where condition ([9]) implies a zero effective dark energy term on the 
brane. In fact, the case W = constant corresponds to the particular case of Randall- 
Sundrum branes. On the other hand, this is the generic class of models one would 
expect from superstring theories, where the size of the volume of compactified extra- 
dimensions is modeled as a scalar field. For example, in low energy Heterotic M-theory 
it is found, after compactifying 6 of the 10 spatial dimensions on a Calabi-Yau 3-fold 
[32] , a superpotential of the form W{(f)) oc e"^, with = 3/2. 

Since in the real world supersymmetry is expected to be broken, it is convenient to 
consider small deviations from the configuration of Eq. ([9]). We do this by introducing 
supersymmetry breaking potentials and f2(0^) at the branes in the following 

way 

U={d^Wy-W\ Xi = W + vu and = -W - V2, (10) 

with |t>i| <^ and \v2\ <^ |iy(0^)|. Potentials and f2(0^) parameterize 

deviations from the BPS condition Q. The precise mechanism by which they are 
generated is out of the scope of the present work. We refer to [HH, [65] for discussions 
on this issue. 

2.2 4-D covariant formulation 

Given the topology M = x 3^/1^2, it is convenient to decompose the coordinate 
system X"^ into {x^, z), where with /i = 1, ■ ■ ■ , 4 covers the foliations and sur- 
faces El and S2, and where z covers the 3^/7^2 orbifold and parameterizes the 4-D 
foliations. With this decomposition, it is customary to write the metric line element 

ds"^ = QABdX^dX^ as 



ds^ = N^dz^ + g^u{dx^ + N^'dz){dx'' + N'^dz). 



(11) 
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Here, N and A^'' are the lapse and shift functions for the extra dimensional coordinate 
z, and g^i, is the induced metric on the 4-D foliations with a ( — h ++) signature. At 
the boundaries we have gj^J = g^u{zi) and g'^fj = g^v{z2). It is possible to show that 
the unit-normal vector n'^ to the foliations has components 

= {-N^/N, 1/N), Ha = (0, A^). (12) 

Additionally, it is useful to define the extrinsic curvature K^^, of the 4-D foliations as 

= ^ [9',u - V^AT, - V.AT^] , (13) 

where ' = dz and covariant derivatives are constructed from the induced metric 
gfj,u in the standard way. Another way of writing the extrinsic curvature is 2NK^i, = 
£Nng^J.u, where is the Lie derivative along the vector field Nn^. 

The present notation allows us to reexpress S'buik of Eq. ([6]) in the following way 

^buik = ^ f dz !d'x^N(R - \K^,K^^ - K'] - h<p'/N)' 

+ {g^'' + N^'N''/N^)d^<pd,<P - 2N-^N'' d^<P 9,0 + f/]) , (14) 

where R is the four-dimensional Ricci scalar constructed from g^^ and K = g^'^K^^, 
is the trace of the extrinsic curvature. Observe that the Gibbons-Hawking boundary 
term Sqr, which appeared originally in S'buik, has been absorbed by the use of metric 
(|TT|) . Let us clarify here that the integration in Eq. (I14p along the fifth- dimension is 
performed on the entire circle S^, instead of just half of it. We should keep in mind, 
however, that degrees of freedom living in different halves of the circle are identified 
through the Z2-symmetry. 



2.3 Dynamics and boundary conditions 

In this section we deduce the equations of motion governing the dynamics of the fields 
living in the bulk and the branes. These equations are obtained by varying the total 
action of the system S'tot with respect to the bulk gravitational fields A^, A^'^, g^j^i, and 0, 
taking special care on the variation of the boundary terms. The brane tensions Ai and 
A2 and matter fields and \l/2 play a decisive role in determining boundary conditions 
on the bulk gravitational fields at Si and S2, respectively. The variation of Stot with 
respect to A^, A^^ and g^i, respectively, gives 



W^(f) - U 
-—(b'd,,c 



(15) 
(16) 
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+3KK^, - 2K^^Kl - ^g,, [{<P'/Nf + g''''d,<pd,<l) + U] + ^d,<pd,<p. (17) 

To write these equations we have adopted Gaussian normal coordinates, which cor- 
respond to the gauge choice A^'^ = (with this gauge one has X^n = d^). Here 
Gfj,y = R^u — \gnuR is the Einstein tensor constructed out from the induced metric 
g^ij. The bulk scalar field equation of motion can be deduced either from the previous 
set of equations (by exploiting energy momentum conservation), or just by varying the 
action Stot with respect to 0. One obtains 

[V^ {Nd,<P) + (07iV)' + K<j>'] = (18) 

Near the branes the variation of the action leads to a set of boundary conditions known 
as the Israel matching conditions [66]. In the present model, they are given by 

K,. - 9,uK = g,, - \m^%,, (19) 
0' = iV5^Ai, (20) 

at the first brane Ei, and 

K,, - g,,K = -h., g,, + ^M^'^T^^, (21) 
0' = -iV9^A2, (22) 

at the second brane S2. In the previous expressions we have defined the 4-D energy- 
momentum tensors T^^ and T^^, describing ^E'l and \E'2 in the conventional way 



2 5Sr 



„9 2 5S' 

and T, 
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'''' V^^g^" 



(23) 



Z2 



where 5*1 and S2 are the terms appearing in the action ^branc of Eq. ([8]) 



2.4 BPS solutions 

Let us, for a moment, assume that the supersymmetry breaking potentials fi(</'^) and 
f2(0^) defined in Section [2.11 and brane matter fields and ^^2 are absent. Then, 
given a superpotential Vr(0), the bulk scalar field potential and brane tensions become 

U = {d^Wy - Ai = W{(f)^), and A2 = -W{(f)'^). (24) 

Under these conditions the system presents an important property which shall be 
exploited heavily during the rest of the paper: The system has a BPS vacuum state 
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consisting of a static bulk background in which branes can be allocated anjrwhere, 
without obstruction. Indeed, suppose that the bulk fields depend only on z, and write 
diiu = ^'^{z)Ti^y, where t]^^ is the Minkowski metric, then one finds that the entire 
system of equations (fT5|) -(fT8l) are solved by functions uj{z) and 0(z) satisfying 

— = -^IV and <p' = NdsW. (25) 

Remarkably, boundary conditions (fT9l) -( j22l) are also given by these two equations. 
Thus, the presence of the branes forces the system to aquire a domain-wall-like vacuum 
background, instead of a fiat 5-D Minkowski background. This property allows us to 
handle the complicated system of equations ( |T5|) -(|T8l) by hnearizing fields about this 
state. This will be considered in detail in Section El 

Notice that the warp factor uj{z) may be solved and expressed as a function of (f){z) 
instead of z 

1 



uj{(f)) = exp 



4 



a 



where = ^ . (26) 



2.4.1 Dilatonic braneworlds 



In the case of dilatonic braneworlds one has W = A e""^, where A is a mass scale 
expected to be of order M5, and a is a dimensionless constant. In this case the relations 
of Eq. (!25l) permit us to solve the background values (p and u in terms of z. Using the 
gauge = 1 for definiteness and assuming A > 0, one obtains 

(f){z) =<f),--\n[l-a^Woz], (27) 
a 

cu{z) = [l-a'WozY^'''\ (28) 

where (pi = (j){0) and Wq = Ae""^'^. Notice the presence of a singularity u = at 
z = l/a'^Wf). Without loss of generality, one may take the position of the first brane 
El at 2; = (since A > 0, this is a positive tension brane). Then, the second brane S2 
can be anywhere between z = and z = 1/a'^Wo. Later on, we shall study the case in 
which E2 is stabilized at the singularity. 

2.5 Effective theory 

To finish, we present the effective theory describing the dynamics for the zeroth-order 
fields from the 4-D point of view. The effective theory is a bi-scalar tensor theory of 
gravity of the form [521 [51] 

S = J d''x^[^R - ^g^'^rX^ad.^b - V 

+Si[^i,Alg,,] + S2[^2,Alg,,l (29) 
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where uja-, with a = 1, 2, are the values of the warp factor uj{z) at the brane positions 
Za and = IGnM^. Observe that Ua can be expressed in terms of 0" (the boundary 
values of 0) by using Eq. fpU|) evaluated at = 0". The elements of the sigma model 
metric 7"* are given by 
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6M5 



1 



where Al 



2Mr,A{ 

Wi 

uj^/B"^, and 



2-, 



7 



22 



6M, 



5W2 



1 + 



W2 



21 



21 



do; 



1^' 



Finally, the effective potential V is found to be 

V{c^\c^') = ^-^[Atv, + Alv,]. 



(30) 



(31) 



(32) 



This effective theory can be deduced either by solving the full set of Eqs. ( IT^ - ffTS|) at 
the linear level [51], or by directly integrating the extra-dimensional coordinate z in the 
action ^ using the moduli-space-approximation approach [H2]. We should mention 
here that Newton's constant, as measured by Cavendish experiments on the positive 
tension brane, is given by Gn = G^A\ = G^uf/B^. 



3 Linearized gravity 

In this section we deduce the equations of motion governing the low energy regime of 
the system -close to the BPS configuration presented in Section [?!^ and consider the 
problem of stabilizing the moduli. Our approach will be to linearize gravity by defining 
a set of perturbation fields about the aforementioned static vacuum configuration. 



3.1 Low energy regime equations 



To start with, assume the existence of background fields (po, Uq and Nq, depending on 
both x'^ and z, and satisfying the following equations 

= -}-NoW{<Po), and 0[, = Nod^Wo- 

The bulk scalar field boundary values are defined to satisfy (pKx) 
4>l{x) = (f)o{x, Z2). The form of the warp factor c^o is already known to us 



(33) 

60 (x, 2:1) and 



ojq{z, x) = exp 



00 



(34) 



10 



where 0o is an arbitrary constant. Now, we would like to study the system perturbed 
about the BPS configuration of Section [231 To this extent, we define the following set 
of variables h^^, ip and as 



N = Noe^, 



(35) 
(36) 
(37) 



where g^^, (j) and N satisfy the equations of motion ffT^ . ffTTj) and ffTSj) . taking into 
account the presence of matter in the branes and the small supersymmetry breaking 
potentials vi and V2. Additionally, ^^i, is defined to depend only on the spacetime 
coordinate x. The functions h^y, ip and ip are linear perturbations satisfying \h^y\ <C 
ujl\gfj_u\, \ip\ <^ 1 00 1 and \ip\ <^ 1. Now, if we insert these definitions back into the 
equations of motion ( fTSl) . (|T71) and ( ITSl) . and neglect second order quantities in h^j^u, ip 
and if) we obtain the required equations of motion for the low energy regime: First, 
Eq. dTHD leads to 



M^(0o) 



(900 



iVo(Xo + X). (38) 



Equation (|T71) leads to 

dNodW 



Ar2 



500 2 900 



(39) 



And finally Eq. ( ITSl) gives 



-A^o 



A^oW^o + 
-7/; - 



2^2 



2— — -0 + ^ttfV^ 



■ft' + 



{Zo + Z). 



(40) 



In the previous equations Uq = t/(0o) and Wq = W{4>o). Notice that the trace h = 
g^'^h^y is taken with respect to g^j^u instead of g^j^u- Equations (j38l) and (!39l) correspond 
to the linearized 5-D Einstein equations, while Eq. (HOl) corresponds to the linearized 
bulk scalar field equation. Notice the appearance of the sums Xq + X, Y^^ + Y^^ and 



Zq + Z aX the right hand side of Eqs. (!38l)-( l40|) . The quantities Xq, and are 



Xo 



(V0o)'-^/? + 8cUo-^ncuo, 



(41) 
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y: 



Gfj,u + - 



-3duUJod^{NoLJo) 



iVr 



whereas X, Y^^^ and Z are 



X = ijQ^-{nh-v"v'^h 



Y 



af3) 



(42) 
(43) 

(44) 



(45) 
(46) 



Operators such as V and □ are constructed out of g^j^i, instead of gfj_y. In writing 
X, Y^y and Z, we have neglected terms involving products between background fields 
spacetime derivatives, such as V ^j,uJo, and perturbation fields spacetime derivatives, such 
as V^/i. This is justified as we shall later consider the stabilization of the background 
fields. Boundary conditions ffTI?l) - fl^^ can also be expressed in terms of linear fields. 
At the brane E^, with a = 1,2, they take the form 



dW 

w{(i)o)iJ + j--ip 

d(po 



and 



dW , d'^W dva dva , d'^Va 



'JO 



^0 



(47) 



(4J 



where signs =F stand for the first and second brane respectively. Background quantities 
like 00 and Nq must be evaluated at 2; = according to the brane. It is also useful to 
recast Eq. ( JT6l) in terms of linear variables 



2iVo 



— (v^V- v^^) 



3 ~ dW ~ 



(49) 



3.2 Gauge freedom 

It is important to keep in mind that Eqs. (|T^ - (|T^ were written in a gauge = 0. 
This gauge is appropriate for studying parallel branes as we are in the present case. In 
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general, given a set of small arbitrary parameters ez{x, z) and e^(a;, z), it can be shown 
that the perturbed theory is invariant under the following set of gauge transformations 



TV' 



h 



(50) 
(51) 
(52) 



The gauge parameter e^{x, z) can be used to eliminate A*"^ from the perturbed theory 
as we have done. The gauge parameter ez{x,z) can be used similarly to redefine (or 
eliminate) ip. Observe that there is a residual gauge freedom to choose e^{x, z) without 
spoiling condition A^"^ 



0. Indeed, if e^{x,z) satisfies 



(53) 



then we can redefine h^y and continue keeping A*"^ = 0. This gauge freedom makes zero 
mode gravity invariant under diffeomorphisms, as it should be. 



3.3 Homogeneous solutions 



Observe that the most general set of solutions h^y, ip and if can be written in the 
following form 



h 



flU 



h^u + h 



1p = 1p + 



and 



if = (p + (f. 



(54) 



Here, fields hf^^, ip and (f are the specific solutions to the system, i.e. those solutions to 
Eqs. ( l38l) -( H0l) and boundary conditions ( H71) -( l48|) including the inhomogeneous terms 
Xq, Y^^ and Zq. On the other hand, fields h^y, ip and (f are homogeneous solutions 
satisfying Eqs. (l38l) - (H0l) but only with X, Y^^y and Z at the right hand side. They also 
satisfy the following linear boundary conditions 



h' 



9f^uh H -— [h^y 



dW 

w^(0o) i^ + i^v 



idVa ^ 



+^Noujlvag^yi! + ^Nou^'^g^yif, 



(55) 



and 



(90^ d(f)o 



^0 



(56) 



at both branes a = 1,2 respectively. Observe that matter fields do not appear in this 
set of boundary conditions. It was shown in [51] that the specific solutions hfj,y, ip 
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and are related to the zeroth-order fields in a special way: They are generated by 
the evolution of the zeroth-order fields ujq{x^ z), No{x, z) and (po^x, z) on the bulk and 
branes and, when integrated, they give rise to the effective theory shown in Section 
12.51 In this article we are utterly interested on the homogeneous solutions h^^, and 
(f. They appear linearly coupled to the matter energy momentum tensor T^j, in the 
brane, which is just what we need to compute corrections to the Newtonian potential 
at short distances (see Section H]). 

3.4 Stabilization of the moduli 

It is clear from the effective theory shown in Section [23] that, in the absence of super- 
symmetry breaking potentials Vi and V2, the scalar fields 0i and 02 are massless. Recall 
that 01 and 02 are the boundary values of the bulk field at the branes (we could 
have equally chosen a combination between the radion and only one of the boundary 
values, say 0i). Solar system tests of gravity provide strong constraints on the con- 
formal couplings Ai and A2 between the moduli and matter fields (recall Section [275|l . 
at the extent of making difficult to reconcile natural values for the parameters of the 
model and observations |l2]. For example, in the case of a dilatonic superpotential 
IV(0) oc e°"^, solar system tests require < 1.5 x 10^^, whereas in 5-D Heterotic 
M-theory one expects = 3/2. For this reason, we consider the stabilization of the 
moduli by introducing boundary supersymmetry breaking terms Vi and V2, implying a 
potential 

V{<p\<p') = ^-[Atv2 + Atv,]. (57) 

To be consistent with low energy phenomenology, we shall further assume that the mod- 
uh are driven by this potential to fixed points such that ^1(0"*^) = t'2(0^) = d^vi{(f)^) = 
<9(/,f2(0^) = 0, implying a zero effective cosmological constant|l| Under these conditions, 
the zero mode fields 0^ and 0^ acquire masses proportional to 9^fi(0^) and 9|f2(0^), 
respectively. On the other hand, the small field ip appears coupled to the branes also 
through terms proportional to d'^Vi{(l)^) and d'^V2{(f)'^) ■ As we shall see in the following, 
the presence of these couplings drives the system to a stable configuration in which 
scalar perturbation fields satisfy (p = h = 0, and only the traceless and divergence-free 
part of h^i, is free to propagate in the bulk. 

To show this, let us start by fixing the gauge -0 as 

■0 = —aip, (58) 

^The pair of conditions vi{(j)^) = and f2(0^) = are not strictly necessary, as present cosmological 
observations indicate the existence of a non-negligible dark energy term in our universe. 
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and define the traceless graviton field 7^jy = hf^^y — ig^^uh. With these considerations in 
mind, the homogeneous equations of motion become 

W 



h! H h 



2dW , 2 



dU 



^=^{Uh-\v^V^l.p), (59) 



{Uh - ^V"V^7, 



a/3 J 



1- 

3 

+2aN'^U^, (60) 



7: 



1 



-aAr2(-^^,n(^-2V^V.<^), (61) 



NW + 



ON dW 



, lN^dW_ 
2~ 



— — - 2a— 



^ = ^Dy.. (62) 



Additionally, boundary conditions acquire the form 

, NW 

NW 

2 

,,,, 92^/ 91^ d'^Va 



(63) 
(64) 
(65) 



at both branes [recall that we are using fi(0^) 



V2[ 



9^i;i(0i) = 9^^72(0') = 0]. 



In the present gauge, Eq. (H9ll becomes [(^^70-/^ — |V^/?,)/ci;^]' = 0. This means 



3^ 

4 



(66) 



where f^{x) is some vector field independent of coordinate z. Notice that the z- 
dependence oc uj"^ of the combination V°"7cr^ — fV^^/i is the same one of the zero-mode 
graviton. This allows us to absorb V^'ja^i — f V^/i in the definition of the zero mode 
graviton uj'^g^u, and set f^{x) = everywhere in the bulk and branes: We achieve 
this by exploiting the remaining gauge (1531) . Then, by evaluating Eq. (1591) at the 
boundaries, one obtains 



uj'aW- 



uj'aW- 



0. 



(67) 
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Thus, unless and are zero, (p and (p' must be null at the boundaries. This 
forces the perturbation field if to stabilize in the entire bulk. Strictly speaking, this 
argument is only vahd for > |^^|. For small values <^ l^^l' '^^^ to 
take into account higher order terms in the expansion of 0, and gfj,^, and the first 
order perturbation (p would not be stabilized at scales of phenomenological interest. 
This is also true for the case in which the squared mass of the (p excitation is larger 
than ~ Wq^^. With (p stabilized, it is easy to check that the graviton trace h 
and divergence V^'y^u are also reduced to zero. The only mode not affected by the 
boundary stabilizing terms Vi and V2 is the traceless and divergence- free tensor ■y^^, 
which is left satisfying the following equation of motion 

- 1^^^- + ^Ar^[2(5.W^)^ - W'h,. = -^57.., (68) 



and boundary conditions 



%u + = 0. (69) 



In the next section we consider solving this equation and show how 7^1, introduces 
modifications to general relativity at short distances. 



4 Newtonian potential 

Now we consider the computation of the Newtonian potential for single brane models, 
in which the second brane S2 is assumed to be stabilized at the bulk singularity a; = 0. 
The tensor 7^^, described by Eq. flU^ has five independent degrees of freedom. From 
the four dimensional point of view, these are just the necessary degrees of freedom to 
describe massive gravity. In fact, as we shall see in the following, there is an infinite 
tower of massive gravitons with masses determined by the boundary conditions at the 
orbifold fixed points. First, notice that Eq. ( l68ll can be further simplified: Assume 
a 4-D Minkowski background, and let 7^1, = e^^'^co^^'^^miz) with = — m^; consider 
also the gauge N = oj. Here $m(^) is the amplitude of a Fourier mode representing a 
graviton state of mass m (for simplicity, we are disregarding tensorial indexes). This 
leads us to the following second order differential equation 

[-dl + v{z)\^^ = m'^^, (70) 

where v{z) is given by 

v{z) = -u^W^ 
^ ' 8 



a 



(71) 
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The boundary conditions for ^rn{z) are now 

C + lujW<!>m = 0, (72) 

at both branes. Notice that ^miz) defines an orthogonal set of fields. Indeed, from 
Eq. (1701) it directly follows 

= [^m{z)Ki^) - KX^)^n{z)] , (73) 

m^' — 

relation which, after integrating and applying boundary conditions (!72|) . gives 



Z2 



dz^rn{z)^n{z) = Snm, (74) 

provided that the $m's are correctly normalized. To solve Eq. (!70l) it is necessary to 
know the precise forms of W and uj as functions of z. They, of course, must be solved 
out from the BPS relations 

u' 1 dW 

^ = - ujW, and (P' = uj^. (75) 
4 ocp 

The first of these two equations gives 



4 



z 







u;'= [1 + ] I Wdz] , (76) 



where we have imposed ci;(0) = 1 and assumed, without loss of generality, that the first 
brane is located at 2; = 0. Notice that there is a rich variety of possibilities for the 
function v{z), depending on the form of the superpotential M^(0). In Section IT2] we 
shall focus our efforts on the simple case of dilatonic braneworlds W{(j)) = A e""^. There 
we find that, depending on the value that a takes, one may have either a continuum 
spectra of massive gravitons, or a discrete tower of states. 



4.1 The potential 

We now compute the short distance effects of bulk gravitons on the Newtonian po- 
tential. For this we consider the ideal case of two point particles of masses mi and 
7712 at rest on the same positive tension brane. To proceed, observe that the traceless 
graviton field 7^,^ comes coupled to the brane matter fields through the term 

Ant = -^i.^T'^-'Siz), (77) 
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where S{z) is the Dirac delta function about z = 0. It is then possible to show [H] that 
the Fourier transformation of the Newtonian potential V{r) describing the interaction 
between two sources with energy momentum tensors T^" and T2'' , is given by 



V{k) 



2M| 



'^1 ^Ul 



(m) rpaf3 



(7J 



where $m(0) are the normalized graviton amplitudes evaluated at z = 0, and Pj^^fB 
the polarization tensor for the graviton mode of mass m [67]. This result comes from 
the Kallen-Lehmann spectral representation of the graviton propagator. In the case 
of point particles of masses mi and m2 at rest one has Tj^'^(k) = mj^Q^Q. This means 
that the only relevant components of the polarization tensor are the 0000 elements 

2/3 for the case of massive 
8]) and Fourier transforming 



(0) 



0000 



1/2 for the case of massless gravitons, and -foooo 
gravitons with m > 0. Putting all this together into Eq. 
back to coordinate space, we obtain 



V{r) 



1 171x1712 



r 



(79) 



m>0 

Observe that the zero mode amplitude satisfies |$oP = M^^^"^ (0) / B'^ , where B"^ was 
defined in Section 12. 5[ This gives the right value for the Newtonian constant G^^ = 
IGnM^B"^ as defined in the effective theory for the zero mode fields. We thus obtain 
the general expression 



V{r) = -G.^ [1 + fir)] 



10) 



where /(r) is the correction to Newton's inverse-square law, defined as 

4, 



/w = T.i$o(o)r^$^i<i',„(o)|v 



(81) 



m>0 



4.2 Newtonian potential for a = constant 

It is possible to compute an exact expression for /(r) in the case a =constant (dilatonic 
braneworlds) . The Randall- Sundrum case is reobtained when a = 0. For concreteness, 
let us take W{(f)) = A e"*^ where A > is some fundamental mass scale. Then, solutions 
to Eqs. ( !75l) are simply 



uj{z) 



W{z) 



'l+'-^W,z 



Wa 



l+'-^W,z 



12) 



^3) 
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where we have defined Wo = A e"*^^ , with the value of cj) at the positive tension 
brane. The form of v{z) is then remarkably simple, depending on the value of a. If a 
is in the range < < 1/4, then 



[1 + kz)^ 



with 



where k = {1 — 4:a^)Wo/4:. Notice that, although the singularity is at a finite proper 
distance from the positive tension brane (as we saw in Section 12.4.11) in this coordinate 
system the singularity is at 2; = +00 (recall that we took N = u). If = 1/4 then 
tu{z) = e~^o^/^ and W{z) = Woe^°^^^, and the potential v{z) is just a constant 



v{z) 



3Wn 



i5) 



Observe that in this case the singularity is also at infinity. Finally, if a is in the range 
1/4 > then 



viz 



with 



.^^(4a^-ir^ 



1 

2' 



^6) 



where /i = (4a^ — l)Wo/4:. Observe that in this case the singularity is at z = l//i. In 
all of these three cases, one has 



l*o(0)p 



1 + 2^2 



Wo, 



(87) 



which means that G^ 



It is interesting to notice that for the cases 



< < 1/4 and = 1/4 there is a continuum of massive gravitons, irrespective of 
the fact that the extra-dimension is actually finite. This is due to the warping of the 
extra-dimension and the presence of the singularity. In the following, we find solutions 
to this system case by case. 



4.2.1 Case = 1/4 

Let us start with the simplest case. Here v = (3VFo/8)^ is just a constant, and solutions 
are given by linear combinations of trigonometric functions. The singularity = is 
aX z = +00, so it is convenient to keep the second brane at a finite position z = Zs to 
impose boundary conditions, and then let Zg —>■ +00. Normalized solutions, satisfying 
appropriate boundary conditions are 



$m(^) 



z^m^ 



Vm'^ — vcos{Xmz) + sm{Xmz) 
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where A, 



\Jvn} — V. The masses are quantized as 

2 



m 



V + 



^9) 



with n = 1, 2, 3, ■ ■ ■. Equation ( ]89ll allows to define the appropriate integration measure 

V dm, 



over the spectra in the limit Zs —>■ +oo. Indeed, one finds Yin ^ ^ I ^/v^ 
with the integration performed between ^/v and +oo. Then, putting it all back together 
into Eq. ( IHTl) . we obtain 



fir) 



32 



dm 



\fm} 







m 



(90) 



4.2.2 Case 1/4 

Here the singularity is also at 2; = +00, so we use the same technique as before, and 
let Zs +00 after imposing boundary conditions. General solutions to ( 1701) are 



^rn{z) = Vl + kziAmJt 



j{l + kz) 



+ B Y 



j{l + kz) 



(91) 



where Ju{x) and Vj,(x) are the usual Bessel functions of order u. Recall that here u 



1(1 -4a2) 



+ |. Boundary conditions at the first brane position give AmJu-i [m/k] + 



BmYv-i [m/k] = 0, which allows to write 



^rn{z) = N^Vl + kz (^J,„i [m/k] Y, + kz) 



771 1 \ 

Y,_^[m/k]J, ^-{l + kz)\) 



(92) 



Boundary conditions at the second brane z = Zg give 

y.-i [m/k] ^ n_i [(1 + kzs)m/k] 
Jy-i [m/k] Jjy^i [(1 + kzs)m/k] ' 

For Zs —>■ +00, this condition implies a quantization of m of the form m„ 
turn, permits us to define the integration measure over the spectra as J2m ^ ^ I ^^n- 
Now the integration is between and +00. The normalization constant Nm of Eq. 
is found to be 



(93) 

— which, in 



N- 





kzs 
Tim 



m 



dz{l + kz) [J^^i[m/k]Y^ —{1 + kz) +Y^_i[m/k]J^ —{1 + kz) 



k 



m 



k 



{jl,[m/k] + Yl,[m/k]). 



(94) 



The second equality comes out in the limit z, 
f{r) by using Eq. (IHTl) 

8 1 - 4^2 ^00 



37r2 1 + 2a2 m Jl_^[m/k] + Y^_^ [m/k] 
20 



00. All of this allows us to compute 

p—mr 

(95) 



where we used the identity Jjy_i[a;]yjy[x] — Ji,[x] = —2lnx. It is of interest to 

check whether Eq. ( I90ll is reobtained out of Eq. ( l95l) by letting 1/4. This is 

indeed the case: It is enough to use the following identity in Eq. flU^ . valid in the limit 
u — >■ +00 (which is equivalent to ^ 1/4) 



TT 



where 9{x — 1) is the unitary step function about x = 1 



(96) 



4.2.3 Case > 1/4 

Here the singularity is at a finite coordinate z = 1/ fj, and u = |(4q;^ 
solutions to Eq. (JTOll are of the form 



1) ^ — |- General 



+ 5„K 



m 



l-fiz) 



(97) 



Boundary conditions at 2; = and z = l/fi require Bm = and J^+i[m//i] = 0. Let 
be the n-th zero of the Bessel function J^+i[x]. Then, quantized graviton masses 



u 



u+l 



are given by 



rrir 



and normalized solutions are easily found to be 



1 — /iZ 



J. 



m 



fiz) 



{91 



(99) 



It should be noticed that solutions with < 0, in principle allowed by Eq. flTUj) in the 
range 5/8 < a^, are discarded by boundary conditions. This means that there are no 
tachionic states in the graviton spectra, as it should be|f| Then, /(r) is simply found 
to be 



fir) 



4 4^2 



31 + 2a 



,2 / V 



(100) 



Interestingly, this corresponds to a tower of massive states contributing Yukawa-like 
interactions to the Newtonian potential, all of them coupled to matter with the same 
strength. Again, we should check that Eq. fl90|) is reobtained out from Eq. fllOOl) in the 



"''It would be interesting to investigate whether this result can be extended to any type of super- 
potential W{(j)) involved in solutions of Eq. ([70)) with boundary conditions (j7^ . 
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limit o? 1/4. This is possible by noticing that in the limit v +00, the following 
relation involving Bessel zeros u^^t^ is satisfied 



n 



V 



y^iK/i^y - 1 + arctan fl/v/(</z^)2 - 1 



2u - 1 
4 ' 



(101) 



This relation allows us to define the integration measure over the graviton spectra 

4 



E 



TrWoUa^ - 1) 

n 

which is enough to recover Eq. (l90i) . 



dm 



— V 



m 



(102) 



4.3 Short distance corrections 

The only length scale available in the present system, apart from the Planck scale, is 
Wq^. Then, we may ask what the leading correction to the Newtonian potential is in 
the regime rWo ^ 1. Since such correction would be the first signature expected from 
these models at short distance tests of gravity, their computation allows us to place 
phenomenological constraints on the values of Wq. 

4.3.1 Case = 1/4 

In the case = 1/4 one finds directly, by using rWo ^ 1 in Eq. (!90l) 



V{r) 



772,1772,2 



(103) 



At present, we know of no constraints on this type of corrections to the Newtonian 
potential. 

4.3.2 Casea^<l/A 

In the case < 1/4 one can expand the Bessel functions in the small argument limit 
to find 



V{r) = -Gi 



727,1727,2 



r 



1 + 7: 



l-4:a'^B[u-l,u- 1] 
3 1 + 2a2 (2fcr)2{'^-i) 



(104) 



where k = v = \( 



1 - 4a2)-i + i, and B[x, y] = T[x]T[y]/T[x + y] is the 



usual beta function. Constraints on k and a for a few values of 2(i/ — 1) = 1, 2, 3, 
appearing in the power law correction r^^*^"^) can be found in ref. [5]. 
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4.3.3 Case 1/4 

Finally, in the case cP' > 1/4, one may just pick up the leading contributing term 
involving the first root Uq~^^. This gives a Yukawa force correction of the form 



V{r) 



17117712 



r 



4 4^2 

1 + 



3 1 + 2a2 



(105) 



It is interesting here to consider the particular case of 5-D Heterotic M-theory, where 



a 



3/2. In this case, the leading correction to the Newtonian potential gives 



V{r) = -Gi 



7nim2 



r 



5 
3^ 



1 + -e-^" 



(106) 



where A ^ ~ 4.45Wo- Current tests of gravity at short distances P provide the con- 
straint A < 50/im. 



4.4 Extra-dimensions in the near future? 

A sensible question regarding this type of model is whether there are any chances 
of observing short distance modifications of general relativity in the near future. To 
explore this, notice that the relevant energy scale at which corrections to the conven- 
tional Newtonian potential become significant is Wq = Ae"'^^ instead of the more 
fundamental mass scale A. Typically one would expect A ~ M5 which has to be above 
TeV scales to agree with particle physics constraints [HHl EHl [TQl [71] . Nevertheless, the 
factor e°"^'^ in front of A leaves open the possibility of bringing A = Wq^ up to micron 
scales, depending on the vacuum expectation value of 0i. 

In the case of 5-D Heterotic M-theory, for instance, one has e°"^'^ = 1/V, where V 
is the volume of the Calabi-Yau 3- fold in units of M5. In order to have an accessible 
scale A ~ 10/xm, it would be required 

V^I-IO-, (107) 

where we assumed A ~ M5. On the other hand, Newton's constant also comes deter- 
mined by a combination of Wq and M5 in the form G^^ = -j^^M^Wq^ . This implies 
~ M|V. Thus, to achieve the estimation of Eq. (ITOTD one requires the following 
values for M5 and V 

M5 ~ 10-^°Mpi, and V^^^ ~ 10^ (108) 

Although, these figures do not arise naturally within string theory, they are in no 
conflict with present phenomenological constraints coming from high energy physics. 



23 



In particular, non-zero Kaluza-Klein modes coming from the compactified volume V 
would have masses of order lO^GeV. On the other hand, if M5 is of the order of the 
grand unification scale Mqut ~ lO^^GeV, then corrections to the Newtonian potential 
would be present at the non-accessible scale A ~ 10~^°/xm. 

5 Conclusions 

Braneworld models provide a powerful framework to address many theoretical problems 
and phenomenological issues, with a high degree of predictive capacity. They allow a 
consistent picture of our four-dimensional world, and yet, they grant the very appeal- 
ing possibility of observing new physical phenomena just beyond currently accessible 
energies. 

In this paper we investigated the gravitational interaction between massive bodies 
on the braneworld within a supersymmetric braneworld scenario characterized by the 
prominent role of a bulk scalar field. By doing so, we have learned that it is possible 
to obtain short distance modifications to general relativity in ways that differ from 
the well known Randall-Sundrum case. These modifications represent a distinctive 
signature for this class of models that can be constrained by current tests of gravity at 
short distances. 

The setup considered here consisted of a fairly general class of supersymmetric brane- 
world models with a bulk scalar field and brane tensions proportional to the super- 
potential W{(j)) of the theory. For W{(j)) 7^ constant, the vacuum state of the theory is 
in general different from the usual AdS profile. Nevertheless, in order to have signifi- 
cant effects at short distances -say, micron scales- different from the Randall-Sundrum 
case, it was necessary to stabilize the bulk scalar field in a way that would not spoil the 
geometry of the extra-dimensional space. To this extent, we considered the inclusion 
of supersymmetry breaking potentials on the branes. After this, the only relevant de- 
grees of freedom on the bulk consisted of a massive spectrum of gravitons, with masses 
determined by boundary conditions on the branes. 

On the phenomenological side, the main results of this article are summarized by 
Eqs. fllOSp . fllOSp and fllOSp . They show the leading contributions to the Newtonian 
potential within dilatonic braneworld scenarios, which is what would be observed if 
the tension of the brane is small enough. Regarding these results, we indicated that 
it is plausible to expect new phenomena at micron scales without necessarily having 
conflicts with current high energy constraints: In the case of dilatonic braneworlds, 
for example, the necessary value of the 5-D fundamental scale was M5 ~ 10^°Mpi. 
Additionally, in the case of 5-D Heterotic M-theory, where is related to the volume of 
small compactified extra-dimensions, we found that Kaluza-Klein modes are expected 
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to be of order ~ lO^GeV. 
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